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The gravitational radiation-reaction force acting on perfect fluids at 3.5 post-Newtonian order 
is cast into a form which is directly applicable to numerical simulations. Extensive use is made 
of metric-coefficient changes induced by functional coordinate transformations, of the continuity 
equation, as well as of the equations of motion. We also present an expression appropriate for 
numerical simulations of the radiation field causing the worked out reaction force. 
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I. INTRODUCTION 



The reaction force acting on isolated perfect fluids due 
to gravitational radiation emission has been expressed 
by means of Burke-Thorne-like potentials up to the 7/2 
post-Newtonian (3.5 PN) order, which corresponds to the 
seventh order in power of the inverse of the speed of light 
c 0. However, this particular form is apparently not 
convenient for computational purposes as it involves up 
to the seventh time derivative of multipole moments. A 
similar problem happens with the 1 PN radiation field 
which causes the 3.5 PN reaction force, since up to four 
time derivatives do apply therein. 

It is well known, e.g. see Q, that the order of time 
derivatives in the 2.5 PN reaction force crucially depends 
on the chosen coordinate system. The influence of gen- 
eral changes of coordinates onto the metric coefficients 
of a many particle system has been worked out in e.g., 
to the 2 PN order. In generalizing this result to 
perfect fluids, combining it with the choice of a suitable 
set of variables, we shall be able to reduce the 3.5 PN 
reaction force to a fifth time derivative object. Further- 
more, applying the continuity equation and the equations 
of motion will result in expressions where only one time 
derivative remains left, which poses no problems for a 
numerical implementation. In this way we generalize a 
previous work by Rezzolla et al. to the case of nonzero 
mass-multipole moments. We are then able to give the 
complete set of equations describing the fluid evolution 
up to 1 PN + 3.5 PN order in a similar form as in a 
former paper by Blanchet, Damour and Schafer [f|. In 
addition, we present the full explicit expression for the 1 
PN radiation field, adapted to numerical simulations. 

Though the formalism we propose limits to the case 
of adiabatic fluids, it may still be used in a large range 
of astrophysical applications. Particularly it provides a 
natural way to generalize the simulations of Oohara and 
Nakamura on the coalescence of binary neutron stars [|J 
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achieved at the 1 PN conservative and 2.5 PN dissipative 
levels, by adding the 3.5 PN contributions to the gravita- 
tional reaction force. Our equations are also appropriate 
to study the effect of the gravitational damping, includ- 
ing the mass quadrupole [7], the mass octupole, as well 
as the current quadrupole (treated separately in paper 
0), on the evolution of fluid modes in rotating compact 
stars in the case where the fluid viscosity is neglected. 
The role played by the bulk viscosity in reducing the 
Chandrasekhar-Friedman-Schutz instability 0, may 
actually be significant, but it can be temporarily ignored, 
regarding the complexity of the problem. Our model of 
matter will eventually be improved in future works. Let 
us point out finally that the adiabaticity condition only 
means the entropy conservation of fluid particles along 
their trajectories, and does not impose any specific rela- 
tion between energy, entropy and pressure. The equation 
of state linking these three variables can thus be freely 
chosen. Suitable choices 0, ^3] permit to model for 
instance rather realistic stellar collapse processes within 
the present formalism. 



II. EFFECT OF A GAUGE TRANSFORMATION 
ON A METRIC 

In order to determine the coordinate system that min- 
imizes the order of time derivatives appearing in the re- 
action force, we need to know the specific effect of a 
general change of coordinates on the functional form of 
the metric. More precisely, we must extend the "contact 
transformation" investigated in paper || for systems of 
point-like particles to the case where the gravitational 
field is generated by a continuous distribution of matter. 
To determine the order at which we need to operate, we 
notice that the parts of the metric we are interested in 
contribute to the 1 PN acceleration, as well as to the 2.5 
and 3.5 PN dissipative dynamics. Products of terms cor- 
responding to the 1 PN and 2.5 PN level in the equations 
of motion cannot be neglected since they affect the 3.5 
PN evolution. Therefore, the gauge transformation must 
be a priori investigated up to the quadratic order. How- 
ever, we shall see that the linear approximation is suf- 
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ficient provided the difference between the old and new 
coordinates is given in terms of the new variables. Still, 
some formulas, like formula (|4.2|l which shows how two 
successive gauge transformations of order 1 PN and 2.5 
+ 3.5 differs from the transformation induced by the sum 
of the associated change of coordinates, can only be ob- 
tained within the quadratic approximation. Now, rather 
than limiting ourselves to the quadratic level, we shall 
give for completeness most of the relations presented in 
this section to arbitrary high orders. The resulting ex- 
pressions will be longer, but more general and not funda- 
mentally more complicated. They may be useful for fu- 
ture works, although not strictly required in the present 
paper. The reader solely interested in the application to 
the actual problem of eliminating time derivatives from 
the 3.5 PN reaction force may skip this section. The 
reader exclusively interested in the final result may go 
directly to the end of section IIVI f equations l|4.31|l ). 

In a given coordinate grid, the metric g^ LV (Greek in- 
dices run from to 3, and Latin indices from 1 to 3) 
is a function of the coordinate x a — (ct — x°,x = 
x l ) and depends functionally on certain matter fields 
X\ , X2, . • • , Xk, ... In the case of barotropic fluids, for in- 
stance, they can be the mass density p and the 4- velocity 
field i// 1 normalized to unity = —1), or the bary- 

onic mass density p* and the linear 3-momentum density 
Mi, or other sets of relevant variables. We have: 

9^ = g^(x a ,X 1 (y a ) 1 X 2 (y a ) 1 X k (y a )) ; 

in short form, g„ v = g^{x a , X A (y a )). In our notation, 
the variables x a contain the entire dependence on the co- 
ordinates, whereas y a are seen as mere "dummy" quan- 
tities. The fields XA(y a ) are in fact themselves function- 
al of the components of the stress energy-tensor and of 
the metric, XA(y a ) = FA(y a ,T af3 , g a p), even if it is not 
indicated explicitly for simplicity. The important point 
is that a coordinate transformation, which is passive by 



essence, affects the function Xa — XA(y a )- 

Let us consider now the (exact) change of coordinates 

x a =x' a + e%,(x'e,X A (yf } )), 

where e x , is a function of x'° and a functional of the fields 
X' A (y^). Such a transformation acts simultaneously on 
the coordinates, on the metric and on the fields X A (y a ). 
As a consequence, the new components of g' are related 
to the old ones by 

g'^(x' a ,X' A (y a )) = g^(x a ,X A (y<*)) 

+ 2g K ,(x^XA(y^)dl ) s x x ,(x^,X' A (y f3 )) 

+ g Xp (x a ,x A (y a ))d' fl e*,(x'e,x A ( y e)) 

xd' v e» xl {x'\X' A W)), (2.1) 

where & x represents the partial derivative with respect 
to x , while the brackets around non- underlined indices 
mean their symmetrization. The fields X' A (y a ) entering 
the argument of read X' A (y a ) = F A (y a , T' a P , g' af3 ). 
Since the dependence on T /Q ^(z 7 ) and g' a p(z^) is purely 
functional, the position variables z 1 appearing there must 
be regarded as dummy, but their other arguments arc 
combinations of tensor components referring to the cur- 
rent coordinate system: T' afi = T' af3 (z 1 , g' s , . . .) and 
Sa/3 ~ 9' a p( zl 1 X' B ) ■ In the post-Newtonian framework, 
the computation of g' is actually restricted to a finite 
order in powers of e", if we assume that the latter quan- 
tity is of order of the inverse of the speed of light c or 
higher, namely e x , = 0(l/c). Therefore, when the func- 
tion g lll ,(x a , XA(y a )) is known, the dependence of g' 
on x' a can be obtained explicitly by expanding the right- 
hand side of equation (|2.1|) in powers of e", = x a — x' a 
within the required precision. Since the order of expan- 
sion can be arbitrarily high, we may symbolically write 



g^(x' a ,X' A (y a )) = g,u(x ,a ,X A (y a )) + ^ a ,X' A (y a ))d' x g iat (^,X A (^)) + 2g x{ ,(x^ 1 X A (y^)dl } e x x ,(^,X' A ( y ^) 

+ 00 , \ 1 ...X k \p 

+ E a ( {k £ ; m + 2) ^wm^^x^)) 

+ e^■■■ A H^^^(y Q ))aA 1 ...A fc .gAp(^^^(/))9;4^,(x'^^(^))). (2.2) 



Note that according to the Einstein summation conven- 
tion, repeated indices are summed over all their possible 
values; e x }' " Afc is a short form for the product e x ) . . . e A , fc , 
and d[ n stands for d[ ... d[ . Again, the infinite sum 

Ai...A fc Ai A fe o ? 

must be regarded as a formal series and does not need 



to converge. It has to be truncated at a given order con- 
sistently with the general approximation scheme. There 
only remains, at this stage, to express the metric compo- 
nents appearing in the second member of equation lj2.2|l 
with the help of the new fields X' A (y a ) rather than the 
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old XA(y a )- It can be most easily performed by using 
the concept of Frechet derivative. Let us consider that 
g^(x a , XA(y a )) is a function of the field variable X\, on 
the afhne space K, having its domain on an afiinc space 
£\ (for instance, the space of smooth bounded real func- 
tion operating on M 4 , embedded with a norm, say | |oo). 
By definition, the Frechet derivative of g^{x a , XA(y a )) 
with respect to Xi is, when it exists, the continuous linear 
form, on the vector space associated to £\, that satisfies 

g^ix^Xito"), X A>1 (y a )) - g^(x a ,X A (y a )) 

= £>x 1 [5^(^,X A (^))].^,X 1 (/)+ (|^,X 1 (/)| 2 oo ) 

for any infinitesimal variation 6* / Xi(y f3 ) = X[(yP) — 
X\(yP) £ £\. Following usual notations, we have inserted 
a dot immediately before the variation 5* l X\ {y 13 ) to in- 
dicate the action of the form resulting from the deriva- 
tion. Denoting by C{£\) the vector space of the continu- 



where (5* , Xb 1 , ■ ■ ■ , 8* , X Bl ) denotes the vector to which 
the /-form Dx Bl ...x B Idnv] applies. It tends toward zero 
as c goes to infinity since X' B and Xb coincide when 
e", = 0. As initially required, all quantities entering the 
above expression depend implicitly on the same variables, 
say x a and X A . In the grid x' a = x a - e%,(x'P ', X' B ), the 
new components of the metric are given, at last, by the 
function g'^{x«*,X> A ) = g^(x' a ,X' A )+6* ex/ g^(x' a ,X' A ). 
In the same way, after a further transformation x' a — 
x" a + e'§, {x" 13 , Xg ), the metric components become 
g% (x" a ,X'X)= g^ (x" a >X'£) + S* L// g^ {x"<* ,X'£). This 
proves the interesting relation: 



x» = x"»+e'^{x" a ,X'i) 
+ e^(x" a +s'% l/ (x"f } ,X^ ) X A fyX'X) 

= x^ + e'^'' a ,X'A)+e^(x' ,a ,h) 



ous linear form on £\, we have Dx 1 [g l _ i v(x a , XA(y a j)] £ 
C(£i). Now, C(£i) can be itself embedded with a struc- 
ture of affine space as the real number set R, so that 
we can define the Frechet derivative of functions / : 
£\ — > C{£\) in a similar way as before. The bilinear 
form Dx B [Dx x [g^v{x a 1 XA{y a ))}\ is said to be the sec- 
ond Frechet derivative of g ^(x 01 , X A {y a )) with respect 
to Xb and Xi, and similarly for higher orders. As the 
derivative of the difference g^ v (x a , X[(y a ), XA>i{y a )) — 
g^u{x a , XA(y a )) with respect to x@ or Xs(yP) is equal to 
the difference of the derivatives, the operator Dx x ...Xa 
(k G N) is evidently symmetric and commutes with the 
space-time derivatives d\ 1 ...x h - We are then in order to 
achieve our goal of removing the dependence of the metric 
on X A (/) by expanding g'^(x a , X' A (y a ) - 5*^X A (y a )) 
around the fields X' A (y a ) after a generalized Taylor the- 
orem ^3|. This yields the fundamental relation: 



(2.3) 



l,k>0 

l+k>l 

■(6* <i X' Bl ,...,5* <t X> Bl ). (2.4) 

The preceding transformation between the initial and fi- 
nal coordinates defines a function e x „ by x^ = x"^ + e x „ . 
Equation (|2.4|) shows that this function is not the mere 
superposition of e^, and e^„ but also involves non-linear 
contributions. We already know that the metric compo- 
nents g'^ obtained after the first coordinate change are 
given by g'^ = g^ v + S*^g^. The ten functions g"^ 
can be similarly deduced from g' ■ g" = g' + 6*, g' 

On the other hand, as the grid {x"^} is related to the 
old coordinate system {x 11 } by x a = x" a + e"„, we have 
g" = g^ u + 5* n g iiv (which can be checked explicitly). 
Finally, the quantities g'^ v = g^ u + S* tl ^ g^ repre- 
sent the components of the metric in the coordinates sys- 
tem {x"' a } such that x a = x"' a + e%,{x'" l3 ,X l Jl{yP)) + 



i=i 



' " ( — i)' 

,g^ = g'nu - 9»v = E Dx Bi ...xb, [g^-iSe^XB!,- ■ ■ ,s* x ,x Bl ) 



...Xb 1 [9fJ>v 

] • (K X , X B!,- ■ .,S* x ,X Bl ) + 2D Xbi ...x Bi [ffA(/*] • K/b,,' 



E 



Ai-.-AfeA/) 



-d\ 1 ...\ k \ P Dx Bl ...x Bl [g^u] ■ (K x , X Bi, ■ ■ ■ ,K m , X Bi) 



, U*+l)(* + 2) 
2 

+ fe , 1 g x'"' kP d\ 1 ...\ kP Dx Bl ...x Bl [g\(jj] ■ (SI^Xb^. ■ ■ ,S*^X Bl )d v )£ x , 
+ e x J- Xk d Xl ... Xk Dx Bl ...x Bl [gx P ] ■ (5lX Bl ,...,5lX Bl )d^,d„e p x }j , 
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e%,(x"' fi ,X%(3/P)). They refer to the same gravitational 
field as g' , or g" . In all what precedes, the new coordi- 
nates are defined by implicit equations for given e x , and 



e'£,, e.g. x' a = x a -e°,(x'P,X' A (yP)). However, they can 
be formally inverted by applying recursively the Taylor 
formula and rearranging the summation, 



n,m>0 



S=l fc 1 + ...+fc 3 =n 1 A 



z?[e A i ••• A *«a ( „ ) , ( „)X 



X (»-D ,(l) ,(l) 

xe 1 - Afc "-i ...a A(2 , A ( 2 ,£ Ai "\9 i(I) ^ A(1)£ «].(fi Bl „..,^i B j. 



In most of applications, the post-Newtonian metric components show the structure of a sum of products of elemen- 
tary factors, each of them being itself a n-tuple spatial integral of a given function /(y, Zi, . . . , z n _i,t, AT^y, ct)) of 
the coordinates and the fields. With the help of the Leibniz rule, the Frechet derivative D x A \g ^w{x a , X B{y a ))\ can 
be written as a sum of terms matching the pattern: 

• K , Xb , 

where F(x a , Xc(y a )) is a function of x a and a functional of Xc(y a ))- After the definition of Dx A , the latter expression 
also reads: 

J rf 3 yi K x , Xa J d 3 y 2 . . . d 3 y n F(x a , X c {y a ))d XA /(yi , y a , . . • , y„, t, X D ( Yl , ct)) . 

From what we conclude that Dx A [g^{x a 1 Xs{y a ))}-5* ,Xa admits an integral representation of the type 
Jd 3 y 5* ^A^gnv/SXji). The so-defined quantity Sg^/SXA is known as the first functional derivative of 
9nv{x a i Xs{y a )) with respect to Xa- It is also possible to introduce the second functional derivative 5 2 g^ v / SXaSXb- 
By definition, 

d 3 Yl d 3 y 2 S* £ ,X B (yi,ct)5* s ,X c (y 2 ,ct) 9 ^ = D Xb x c [g^(x a , X A {y a ))] ■ (5* e ,X B (y%6* E ,I c (f)). 

I 



F(x a ,X c (y a ))D XA 



d 3 yi ■ ■■d 3 y n /(yi,y 2 , ■ • ■ , y n , t, A D (yi, ct)) 



The case where the non-linear terms occurring in (|2.3() 
do not contribute to 5* ^g^v (or at least to the order 
of gp,v we are interested in) is of particular importance. 
This happens when all products involving two factors 
equal to 5* l Xb or e x i as well as a space or time deriva- 
tive of the metric dg can be neglected. Hence, the n 
PN transformation law is linear whenever C(l/c™ +1 ) = 
0(e 2 x ,dg)+0(e x >5* x ,Xdg)+0((5* x ,X) 2 dg), and equation 
(|2.3|) reduces then to: 

K a ,9iiy = ~Dx B [g^u] ■ K X , X B 

+ e x x ,dxg^ + 2g x ^d y) e x xl . (2.5) 

As a check, let us specify relation (|2.5|l to the impor- 
tant case where the matter variables are (i) the baryonic 
density of the fluid p* = ^—gpu° (with g = det g^v), (ii) 
the fluid coordinate 3- velocity field v l = cu l /u°, and (Hi) 
the entropy per unit mass s, i.e. g^ = g t _ LU (x a , p*, v p , s). 

By using the transformation properties of p (scalar), 
(vector) and \J—g (density), we find at the linear order 



after a straightforward calculation: 

if e°, 

<Tp* = di(p*e l x ,) <9i£°</9* + —d t p* 

c c 

= di(p*e x , -p* — el,), 
c 

,5V = -d t (e x ,v L ) + e x ,d 3 v l - v J d 3 e x , 
- dts x , + ^-9j£> 

= |(-4< + 1 4) - (-4' + £ -fv J )d 3 v* , 

5*s = —d t s + e x ,d l s , 
c 

having made use of the continuity equation to obtain the 
second expression of 5* p* (see section llV)) . The variation 
S*v l is put into a rather compact form by resorting to the 
total derivative d/dt — d t + v l di. By virtue of equation 
(|2.5|) , the new components of the metric read 
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9'nv = 9^u-D Pt [g llu }-5*p !f -D vq [g tlu }.8*v q ~D s [g t , u }-5*s 

+ 4^A^+25A( M ^)4- (2-6) 

If the source is made of one point particle of mass m, then 
p*{y, t) — mS(y — z(i)), and the metric can be viewed as 
a functional of z and v, 

g^(x a ,z(t),v(t)) = — d 3 y ffMt ,(x Q ,y,v(y,i)) p*(y,t) . 

Tfl J 

The contribution of the functional derivatives is 
- m(D Pr [g^] ■ 8*p* + D vq [g^] ■ 8*v q ) 

= -J d3 y {9^{x a ,yMy,t)W p* 

+ D Vq [g^(x a ,yMy,t))} ■ S*v q P*} 
d 3 y {(D y ,[g^{x a ,yMy,m 

+ diV j D vi \g liV (x a ,y,v(y,i))])} • {p^ x , - p*^-4) 
P*D Vi [g^ a ,y, v(y, t))]} ■ ( j t ( - 4 + ~4) 

-(-4 + ^)^) 

d 3 y {(-4, + -e° x ,)D r [g^(x a ,yMy,t))] 



+ + -4 

at c 

x lM^(x a ,y,v(y,t))]}<5(y - z (i)) , 
which shows that 

Dp, [g^ v ] ■ S*p* + D v „ [g^u] ■ S*v q 

= D z i [g^] ■ 5*z l + D vi [g^] ■ ^S*z l , 

where S*z l = -ej,,(z(t), t) + v l (t)e x ,{z(t) 7 t)/c. This is 
in perfect agreement with the linearized boost trans- 
formation formulas |jjj]: ct' — ct — e x ,(z(t),t); z n (t') — 
z*(t)-sl,(z(t),t). 



III. ELIMINATION OF HIGHEST TIME 
DERIVATIVES IN THE METRIC OF A 
GENERAL PERFECT FLUID 

An analytic expression for that part of the metric that 
contributes to the 3.5 PN reaction force arising in an 
isolated system due to gravitational wave emission is 
available in a Burke-Thorne-like gauge which also fulfills 
the harmonicity conditions at 1 PN order for gen- 
eral smooth, spatially compact distributions of matter. 
However, it is not suited for numerical calculations for it 
depends on seventh order time derivatives of certain mul- 
tipole moments. A possible cure consists in performing 
a change of coordinates such that terms of 6* l g M „ cancel 



the highest time derivatives in the relevant part of g^y 
(the functional dependence will not be specified anymore 
henceforth). In order to be able to apply the formula 
(|2.5|) . we first have to make choice of the fields Xa- It is 
convenient to take the mass density a — (T 0Q + T tl )/c 2 , 
the current density at — T 0l /c, and the stress = T lJ , 
following paper All three quantities have compact 
spatial support and are Newtonian at the leading order 
for weakly stressed systems. Moreover, they lead to a 
simple expression for the metric, based on a parametriza- 
tion by means of certain potentials of which they are the 
sources. At the 1 PN order, these are the Poisson inte- 
grals: 



U{*,t) = G(-^- a(y, 

J l x -y| 



Ui(x,t) = G 



d 3 y 
x-y| 



o-i(y^) , 



(3.1) 



as well as the super potential x( x j t) = G jd 3 y |x — 
y|er(y,i). They represent the Poisson solutions of the 
equations AU = didiU = -AirGa, A{7, = —AirGdi, 
Ax = 2U respectively. We shall write U = A~ 1 (-4:irG<j) 
or x — A~ 2 (— 8irGa). It is worth noticing that U re- 
duces to the Newtonian potential when c goes to infinity. 
The terms of the metric that contribute to the 1/c 5 and 
1/c 7 part of the equations of motion take also a partic- 
ularly simple form. They indeed depend on only two 
"reaction" potentials JJ rcac and [/J eac involving the mass 
quadrupole hj(t), the mass octupole Iijk(t), and the cur- 
rent quadrupole Jjj (t) as defined in the Blanchet-Damour 
formalism (see [l4| and references therein): 



C/ rcac (x,t) 



5c 5 



189 yfeV ' 



(3.2a) 



1 



21 



~Jjk r(6). 



the symbol x 13 , or equivalently x <y fe -- > ) being a short 
form for the symmetric trace-free (STF) part of x %3 = 
x l x J x h . . ., and I^j(t), an alternative notation for the 
nth time derivative d n Iij (t) / dt n . These potentials differ 
from V^ rcac and y. rcac used in respectively by terms 
of order 1/c 8 and 1/c 6 which do not play any role at 
our level [j]. We give for completeness the multipolc 
moments entering (7 reac and JJ^ eac as functionals of a 
and av: 



hi = d 3 y[ f 3 a 



14c 2 



l —y kk f 3 d 2 t a - 



20 
21c 2 



(3.3a) 
(3.3b) 
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Jij = d 3 y e km<l y 3>k a,, 



(3.3c) 



with £kmi holding for the Euclidean Levi-Civita tensor 
(in three dimensions). Having the potentials U, U%, x, 
jyrcac ^ an( j jjreac ^ Q our disposal, we can write the metric 

components as 

.900 = - 1 + \{U + U rcac ) + \ [d 2 X - 2C/ 2 - 4f/C/ rcac ] 



+ (ffOo)(6+8) + O (^Jo^J > 



(3.4a) 



.9o» = - ^(Ui + t/r c ) + Gkm)(b+7) + O f i ) , (3.4b) 



3ij 



= 5i 



1 + ^((7 + [/ roac ) 



+ 0»i)(4+6) +C ( ^ 



(3.4c) 



where the indices between brackets indicate the order of 
the term they refer to, so that <?oo(6+8) denotes the full 
set of 1/c 6 and 1/c 8 terms entering goo- 

It should be stressed that the 1/c 2 order in the cor- 
rections to the flat metric entering the right-hand side 
of equation (|3.4|) is not the post-Newtonian order, the 
latter being rather defined as the level of contributions 
to the equations of motion. Since (300) (2) is responsible 
for the Newtonian force, the terms (300 ) (4) j (3oi)(3)i ano - 
(<?ij)(2) that come just after in the power expansion are 
evidently post-Newtonian. When the matter variables 
Xa are chosen in a way they have a non-zero Newto- 
nian limit Xa(o) = Xa + C(l/c 2 ) and derive from the 
trajectory of the fluid in the configuration space inde- 
pendently from the metric, the 1/c 5 and 1/c terms of 
the equations of motion identify with the reaction force 
in the post-Newtonian scheme. Other set of variables of 
the type Xa — Xa(o) + 0(l/c 2 ) with Xa/o) 7^ may be 
equivalently used provided all expansions are performed 
consistently including the Newtonian, 1 PN, 2.5 PN, and 
3.5 PN terms. The coupling between space-time curva- 
ture and radiation being a 4 PN effect, the even orders are 
purely conservative. It is thus convenient to split g^ LV into 
an odd (<7 M „) dd and an even part (<7 M „) C vcn, named after 
the parity of the terms they generate in the equations 
of motion expressed with the help of the variables that 
have been adopted. As an example, (300) (7); (3(h)(6) or 
(ff»i)(5) will belong to the odd part, while (ffoo)(2)> (3o*)(3) 
or (3ij)(2) will enter the even one. In our problem, the 
odd (resp. even) part includes the 2.5 and 3.5 PN (resp. 
Newtonian and post-Newtonian) corrections. The even 
part corrections that are beyond the 1 PN approximation 
do not play any role and can be let unspecified. 

We propose now to conduct a coordinate transfor- 
mation intended to reduce, as far as possible, the or- 
der of time derivatives appearing in the original met- 
ric components. This can be achieved by eliminating 
the potentials U^ c and U™£f containing the highest or- 
ders of derivation. In addition, we require that the new 
gauge identifies with the standard post-Newtonian one, 



di9' 0i -kdo9'u = 0(l/c 5 ), ^-^04-2oo) = 0(i/c 4 ), 
which is itself equivalent to the Arnowitt-Deser-Misncr 
(ADM) gauge [lj| up to the 1 PN order. As the original 1 
PN gravitational field l|3.4|l satisfies the harmonicity con- 
ditions in the near zone (i.e. the region of space including 
the system in which the post-Newtonian approximation 
is valid) modulo 1.5 PN corrections, the functions e x , cor- 
responds at the leading order to the difference between 
harmonic and ADM coordinates. After paper we 
have e° x , = C(l/c 3 ) and s l x , =0(l/c 4 ), which shows that 



non-linear terms of the kind s x } X2 "d\ 1 a 2 ...A.. [g^v] ■ ■ ■ = 
C(l/c 10 ) drop out in equation (|2.3|) . The only non- 
linearities that may remain are those generated by the 
functional variations 5* , a, 5* t Ui, and 8* ,<7ij. To know 
whether they can be neglected or not, we have to de- 
termine the post-Newtonian order of the three latter 
quantities. The linear approximation is certainly suffi- 
cient for this purpose. We compute o-{y' a , g' a g (z 7 , X' B )), 
<TiW a X' B )), and a l3 (y' a ,g' a0 (z\X' B )), using 

the gauge transformation law T 11 " (y' a , g' a p(z 7 , X' B )) = 

T^{y a iga p{z\X B )) - 2d^T^(y a iga p{z\X B )) + 
0(e 2 ) for the stress-energy tensor. After straightforward 
calculations, we get 

5* a = s'lidkcr — -adts®, + -e°,d t o- — -aidiE ,,, 

e m i x C C C 

2 2 

^o-id t e x , ^o- lk dkS x , 

c c 

+ -§<7 fefe di4 + 0(s 2 ) + 0((Sl,g) 2 ) , (3.5a) 
5* ,<7i = — odtE x , — <TjdjE x , + e^djo-j — -aid t e x , 

x ' ' ' ' C 

+ -£ x ,d t Oi - -cradiel, 



1 



+ ^dtel, + 0(e 2 ) + 0{{5l,gY) , (3.5b) 

J) 



St .o-j 



2a(idte 3 J, - 2d k si,a i)k + e x ,d k a 



j)k "I" £ x >°keij 



1 



+ -e° x ,d t a l0 + 0(e 2 ) + 0((S*g) 2 ) , (3.5c) 



which means notably that 6* Xa — C(l/c 4 ) = 
0(e x ,,e x ,/c). Consequently, there cannot be nonlinear 
contributions in 5* ,3^,,, and equation l|2.5|l yields 



(b~* ,3oo)odd 



* fc a;'(6+8) 



:Ud t el, 



x'(6) 



C v>"i c° 

2 2 
+ -(50i)(3)9t4'(5) + -^ £ x'(6) d t U 



+ £ l' (5) ^i(ffOo) (4) + -^ £ x'(5+7) d i U 

2 f 3U 

-V 2 J d " y Jx- A {5: - XA) ^ 



d 3 y 5 (.900)(4) fx * 



5X A 



(S*X A ) (5) , (3.6a) 



(S*,g i)o 



dd 



7 



(5) 



d 3 y 



(K x ,9ij)odd = 2d {l e J x ] (5+7) + ^%4'( 5 ) 

^P^jd 3 y-^-(SlX A ) {5) . (3.6c) 

On the other hand, relations l|3.3|l and i|3.4fl tell us that 
the time derivatives of highest order in g' , namely d® 
and dj, come (i) from the first correction U^ c of the re- 
action potential {7 reac that parametrizes the 3.5 PN part 
of g' 00 , and (ii) from the leading approximation of the 
reaction potential that contributes to the 3.5 PN part of 
9oi> U- ca - c = We shall restrict ourselves to the case 

where the new metric can be identified with a function of 
certain elementary potentials Pi, P2, ■ ■ ■ These potentials 
are supposed to be spatial integrals whose sources are a 
priori products of (i) densities such as cr, <7j, <7y, or their 
derivatives evaluated at the location of the integration 
variable y at time t, (ii) some factors depending on the 
field point x, and possibly (Hi) instantaneous potentials 
of the same kind. Typical examples are U, Ui, or U^ c . 
That the latter quantity is indeed of the required type, 
this can be seen by letting the factor x y and the time 
derivatives go under the integration symbol in the first 
term of relation (|3.2a|) . Following this hypothesis, we 
shall search for e x , as a function of U, Ui, or some other 
potentials P\,P2 ■ ■ ■ Moreover, our choice must be such 
that at least U^ c simplifies with terms coming from the 
gauge transformation l|3.6|) . 

To be more explicit, we come back to equations <|3.6all . 
Together with equations l|3.4|) . they imply 

2 2 

(ffoo)(7+9) = (^00) (9) - -^* £ ^'(6+8) + ^a £ x'(5+7)® 



%0i)(3) 

sx A 



(S*X A ) (5) , (3.6b) 



J) 



_ (6+8) ~- 

(3.7a) 



(ff0i)(6+8) - (^0t)(8) ~ di £ °x'(6+8) 



1 4 

■ -o t e x , (5+7) - -^u l{5) , 



(3.7b) 



2 

(<4)(5+7) = (M(7) + 2%4'(5 + 7) + 



■ij jjreac 



(3.7c) 



For simplicity, the 3.5 PN terms that are functionals 
of e°, (6) (Pi, P 2 , . . .) or e^( 5 )(Pi, Pi,---) and functions 
of potentials involving second order time derivatives, at 
most, were gathered in the matrix h^ v . At the leading 
order, equation l|3.7a|l reduces to 



(So 



00 ) (7) 



--dial, 



x'(6) 



2 



2°*' (5) 



diU 



2 [ o 5U 2 



reac 

(5) • 



As long as U^ c appears explicitly in the gravitational 
field, we shall have to evaluate the fifth order time deriva- 
tive I^p (t) in numerical simulations. We can avoid this 
by imposing that the reaction potential U^° cancels the 
contribution to (<7oo)(7) that is linear in e®,^ or s z x r( 5 \, 
which amounts to demanding that e x , , g j (resp. e^,/^ 5 j) de- 
composes into terms containing time derivatives of fourth 
order at most, plus terms Sx 1 ^)^ ( res P- ^'(m** ^ or 



(part II) 



) satisfying 



2 j £ o (P art ) 
c 



~ d t £ x'(6) 



1 £ 3 



(part I) 
'(5) 



- SU r k (part H) 

~ 72 i d y "^^'(5) 



While cither 2ei, 



(part I) 



Je (P art n ) • 
--x'(5) 

rreac I Jl 



x'(5) 



diU/c 2 or -2/d 3 y (SU/5a) 



■:'].,,-, 'dk<j/c 2 contribute to the cancellation of 



2[/| cac /c 2 , the quantities e l x , v ,^^ or £ x u^ are, by 



contrast, explicit functions of U^ c . If £^'(5) f or 
^'(sT*^ 7^ 0, the fifth derivative l\^\t) still occurs in 
the metric by virtue of the relation g[j — 2d^e x ,^ + 
0(l/c 7 ). This eventuality being rejected, 2U^ c /c 2 must 
simplify with — 2d t £°/lgf It ' ) / c exclusively, i.e. 



'(6) 

2 a c (P art ) , 2 rrrcac _ n 

Therefore, we may construct £ x it 6 \ as the sum of 



U, 



reac( — 1) 
(5) 



/c = -Gx ij I^(t)/(5c 6 ) (antiderivative of 
^(5) ac / c that vanishes when = 0), plus some ar- 

litra 

of order higher than four: 



bitrary function e x , |gx b ' that does not involve derivatives 



-0 

'-x'(6) 



1 



r reac( — 1) 
(5) 



( arb ) 
~-x>(6) 



(3.8) 



As a matter of fact, the choice of £2'(6) together with 

that of £ x n 5 \ entirely determines the maximal order 

n d((g'oo) (7)) of time derivations in (<7oo)(7)- Next, the ex- 
pression for the odd Oi components of the new metric 
(l3~7bT) at order 2.5 PN, 



'^ £ x'(6) + c ^ t£ x'(5) 



(floi)(6) - - 

shows that n d ((g' 0i )(e) 

-0 (P a 
-x'(6) 



= Mdi £ °x>(6)) 

die®,^ 1 ^ cancels one of the terms entering d t e x , (r) /c 



4, unless 

(5)1 



8 



In the latter case, (<7oi)(6) reduces to — cV^,|gj b ' ) + 

^t £ a;'(5) b V c J where eJ;'(5) ^ s an arbitrary function of el- 
ementary potentials. We have then 



(6) J 



r / (arb), 

max[n d (e a ., (6 - ) J 



, j (arb),, 

nd{s X '( 5 ) )J 



Since n d ((g' 0l )(6)) as well as n<j((fl l oo)(7)) are adjustable, 
we shall succeed in further lowering the order of temporal 
derivatives if nd((g'ij) (5)) < 3. To know whether this 

inequality is fulfilled, we solve the equation —die x ,^ Tt ^ + 



d t el 



(part) 



(part) 



(arb) 



~'x'{5) / c - for e^,, (5) _ c x ,^ °x'(5) 
solution that vanishes when iy (t) = reads 



j (part) _ _ o (part)(-l) 



(-2) 



The 



(3.9) 



which leads to odd metric space components of the 
form (^.)(5) = - 2G 4 3) W /( 5c5 )' with 

«rf((5i 3 )(5)) = 3. No function of eg,[g) b) or 4'(5) nav " 
ing a maximal time derivative order smaller than 3 are 
capable to rule out the second term, so that we may take 

(arb) _ a (arb) _ 
£ x'(6) ~ t x'(5) - U - 

Having specified the gauge vector at the 2.5 PN order, 
we can proceed to the next approximation, and deter- 
mine both £^/( g ) and £^.'(7) by using the same method 
as employed above. The matrix h^v, depending on 
none of the unknown functions, is entirely fixed at this 
stage. The highest order temporal derivations again 
arise in the reaction potential. At this level, we have 
M^Jodd) = n d (U^ c ) = 7. As before, we intend to 
make this number decrease by discarding the term U^ c 
from the metric Ij3.7|l . Its complete or partial absorp- 
tion into £ l x i, 7 \ would result in a explicit dependence: 



(Uf'S? c ), which would entail the appear- 



-x' (7)^(7) 

ance of time derivatives of seventh order in the spatial 
part 13.7cfl of g'^ u . Considering that (5^) (7) — 2d^e x ]^ 
contains fifth order derivatives at most, these high order 
time derivatives could not be erased by means of other 
contributions. We would end up to the same value of 
M(g^)odd) as before: n d ((g'^) odd ) = n d ({g' l0 ) (7) ) = 7. 
Because it is not acceptable, we must rather incorporate 
U*^ c into e x , , g N , and set 



n d (4'S?; b) )<5. (3.11) 

The space components of the new 3.5 PN metric are ob- 
tained by inserting the expression above into the right- 
hand side of equation (|3.7c[l . This yields a maximal order 
of temporal derivations equal to n d (e x ,^) = 5, whatever 

functions £°'[g) b ^ and £^.,|~ we have used. We have still 
the freedom to specify the "arbitrary" part of e x i as most 
convenient for us. 

In summary, it is possible to reduce the number 
n d{g'fj,y) from 7 to 5 by performing a suitable (linear) 
gauge transformation. Equations (|3.8|) . i|3.9|) . I|3.10fl . and 
(|3.11|) provide a possible choice for the gauge vector: 





a;' (6+8) 



1 <-\ Tr reac(— 2) 

-d t U,, 



(5+7) 



-a' (5+7) 



diU, 



z(-2) 



(5+7) 



j ri cac 
2 U i(5) 



(3.12a) 



with 



5c 5 ' 

1a«4 5) (*) 



189 ^ k() 



and 



U 



e(-l) 



G 



21 - ?fe v ' 



0I ? 



(3.12b) 



modulo the "arbitrary" part. There exist other inter- 
esting alternatives. For instance, we could let the origi- 
nal metric unchanged at the 2.5 PN level while keeping 
the requirement «d((ff^ l/ )odd) < 5 since neither (<7oo)(7), 
(<7oi)(6)) or G?ii)(5) involve time derivatives of higher or- 
der. The coordinate system defined by means of formulas 
(|3.12l) present the advantage to coincide with the ADM 
one up to the 2.5 PN approximation, as we shall see in 
the next section. 



IV. EXPRESSION OF THE ODD METRIC AT 
THE 3.5 PN ORDER 



'a:' (8) 



-u, 



reac( — 1) 
(7) 



■ £ 



( arb ) 
x'(8) 



n d (e x ,% h) )<6, (3.10) 



n d{{g' f _ lu )odd) going automatically down to 6. The high- 
est order time derivatives actually come from the terms 
-<9i£°, (8) and -4C/[ ( c 5 a ) c /c 3 in (ffo;)(8), but they can be 

eliminated if a part of — die x ,^° + d t e x ,^ exactly can- 
cel them. This happens for 



For computing the reaction force, we must now final- 
ize the form of the gravitational field by fixing the two 
arbitrary functions £°,|g! | b ' > and e x }^\ On one hand, 
the gauge x 1 ^ = x^ — e x , is constructed in such a way 
that g'^ v is precisely the ADM metric (c/^adm at the 
1 PN order: g' = (<7 M! ,)adm, plus a 2 PN deviation. 
On the other hand, (<?oo)(9) depends on the 1 PN compo- 
nents of the metric l|3.4[l in the Burke-Thorne-like gauge 
introduced in paper Now, the non-spatial terms of 
(<?Aiy) adm differ from those of g^ v at this order. We have: 



--x'(J) 



a rrrcac(-2) 



(5) 



= (-1) 



_i (arb) 
~-x<(7) 



(goo) 



ADM (4) 



-u 2 



9 



G?cm)adm(3) = - ^Ai = ~ \ 4Ui + \didtxj ■ (4.1) 

In order to homogenize the expression of the 3.5 PN met- 
ric, we shall choose the "arbitrary" part of e x i so as to 
discard the contribution of (<7oo)(4) (resp. (<?oi)(3)) ancl 
replace it by that of (£oo)adm(4) (resp. (#(m)adm (3)) in 
(^*5^)odd- This is achieved by taking 

q (arb) ,j o r~\ I 1 I £* V \ 

£ x>(&) = £ x'(5) C>i£ x'(3) — ^bI £ x'(3)J ■ [O^XbJ 



dd / (5) 



i ( arb ) 
-x'(7) 



— £ x'(6) u t x'(3) "T" £ x'(5) a i fc x'(3) 
= = e °'(6)^0£x'(4) + 4'(5)^/4'(4) 

-D XB [4, (4) ].^x fl + oCi) , 



with e odd = ( £ 2' (6+8) >4' (5+7))- This statement can be 
checked either directly or by use of formula l|2.4[) of sec- 
tion [nj noticing that the final value of e x >, which reads 



£ x' ~ £ cvcn + £ odd + £ odd^A£even ~~ ^X B [ £ even] ' %B 

'odd 

(4.2) 

if we pose e£ ven = (4'(3)' e x'(4))> is exactly the one we 
would get by applying successively the two coordinate 
transformations x"^ = x^ — e^ vcn (x" a , X A ) and x'** = 



The resulting metric g' is 
9oo 

cr 

2 A a _i 

'(5) 



x'(6) 



^Aidtel, (5) + -4 (6) 9 t [/+ -4, ^(-2^ 
2 /"., SU 



+ c 2 £ x'(5+7) 



d t U-- d 3 y 



sx A 



1 



1 



d 3 y 



c 8 V c 



J(-2£/ 2 
1 

To 



1 



.(s:x A ) {5) + -(...) 



(4.3a) 



/ 1 2 q 2 

? 0i = _ + ~^2 U9 i £ x>(6) + ~^3~ Ud t £ x'(5) 

L 4 r).cJ Lc-? r) 4 

c 3^ aie a;'(5) c 3 £ z' (5) U 3 ^ 

1 /\, SAi ... _. , 1 



(C,^a)(5) + -(•■■) 



+i(-o+ofi), 



7^ 



(4.3b) 



"* rrrcac 
2 U (5) 



' 2 <% £ x'(5+7) 



2 Hi f,3 <5Z7 ..... , 1 



-(^,^) (5) + -(...) 



(4.3c) 



where £z'(g +8 ) and £ % x i/ 5+ ^ are defined in equations 
(|3.12l) . We remind that (#oi)(6) = by construction. 
We may now calculate the functional derivatives ex- 
plicitly. Considering for example the potential U, its 
source contains none of the densities cr, or ct^ , and thus, 
fd 3 y (SU/Sai)S*ai = fd 3 y(SU/Saij) Pay = 0. The 
term jd 3 y (SU/Sa) S*a is deduced from the evaluation 
of the difference U' — U: 



U'-U : 



= G 



G [ a ' (y, t) - G f a(y, t) 

J |x-y| J |x-y| 

''"• V <5V(y,i)+0(ra] 2 ), 



x- y 



hence 



"" ; y ^(C x , CT )(5+7) 

C' / T— ^-T (<5 £ »(5 + 7) 



(4.4a) 



G 



x- y| 
_^y_ 
x-y| 



'a:'(5+7) afc 



5fc(7 <7<9 t £°, 



x'(6) 



^'(6) 



x'(5) 



x'(5) 



(y,*). 

(4.4b) 



The other non-zero functional derivatives present in the 
metric variation i|3.6f) are computed following the same 
procedure. 



yftfli |d 3 y |x-y| [4 (5 ) 5 ^](y>*)> ( 4 - 4c ) 



rf3 y ^(C x ,^)(5) 



= 4G 



d 3 y 
x-y| 



- CT ^4'(5) -^ 9 J'4'(5) 



£ i>(5) d 3 a i 



(4.4d) 



By inserting formulas l|4.4|l into equation (|4.3|l , we get 

an expression for the metric already well adapted to our 
purpose. Notably, all potentials elementary potentials 
but x an d Gjd 3 y |x — y\{8x/5a)8* t a satisfy Poisson- 
type equations. Nonetheless, in order to make easier the 
comparison with paper we shall present an alterna- 
tive form where the functional derivatives are re-written 
as combinations of Poisson-type integrals on the sources 
Xi = a and Xi — a{ exclusively. It can be obtained from 
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an appropriate handling of the integrands in the poten- 
tials (|4.4|l . Originally, they are indeed made of sums of 
elementary pieces of the kind 



_ll»2...1j_a IrrH-l 



y H - 1 )(x 1 ' — y H ) + permutations] 



+ . . . + (-i)V 1 - y ll )(* 12 - v 12 ) ■ ■ ■ - v H ) 



|x- y 



2p-l 



dX A (y,t)de 

a 

2p— l^»i*a...i 



(Id 



l dXA(y,t) x function of time, 



for A = 1, 2 and p, I 6 N. Now, it is straightforward to 
show by recurrence that 



so that every piece decomposes into a sum of terms such 
as ^'ija-J*^! -y^)(x i2 -y l2 )...(x 1 ' -y^^x-y^P- 1 . 
Each of them are next transformed by means of the idcn- 
(x 11 — y H ) + permutations] tity 



[1/2] 

(^-^)...(^-y*0|x-yr=]T 



(a + 21- 2m) ... (a + 2) 



^2m+i--i! I x _ y| a+2 ' 2m + permutations 



valid for a > — 2, where the latter relation results from the formula 



[1/2] 



SU...<, |x - y| a = ]T a ( a - 2 ) ■ ■ ■ ( a - 21 + 2m + 2 )l x - y\ a - 2l+2m [s 1 ^ . . . 

x (x l2m+1 — y l2m +!) . . . (a; 1 ' — + permutations] , 



m=0 



r 



proved by simple recurrence. The square bracket [1/2] 
denotes here the integer part of 1/2. We thus see that 
the initial potentials can be split into a sum of terms 
matching the pattern: 



d 3 y a**- J *fiU...i, |x - y| 2 f- 1 9X A (y ) t) 



ll2 ... tl Jd 3 y Ix-yl*- 1 ^, 



*)• 



Under this form, they only depend on the Poisson- 
like integrals P q = Gjd 3 y |x - yf'-^y, t), P lq = 
Gfd 3 y Ix-y^-V^y,*) as well as P ijq = Gfd 3 y |x - 
y| 29_1 0Xj (y, t). We shall need in particular P = U, P\ = 
X, Pw = Ui as before, but also Xi — Pn an d Uij — Pyo- 
As an example, let us apply the latter treatment to the 
simplest functional derivative Gjii 3 y/|x — y| e k l( ^.d k a. 
We have: 



d 3 y 



= G 



-— / (3) a) 



(a; 1 - y l ) 



dja 



d 



d 



i dy ^~ yl dy- a 
d 3 y da 



y| dy 3 



After integrating by part and observing that d\x — 
y\ a /dy l = — d\x — y\ a /dx l for any real number a, 
we arrive at the equality: Gfd 3 y e k x , ^d k a /\x — y| = 

-2Gl\f \t){-dijX +x l d j U]/(5c 5 ). Processing all the ele- 
mentary potentials in a similar way leads to the impor- 
tant result 



( I \ 4G 7-(3)a 

(,5oo)odd - ~-^f 1 H "MX 



+ ^(8d k xi + ^d H d t P 2 - x k d l d tX ) 
I s , II ... 29 



+ ^4H~x k x l U - ^r 2 d klX - ^x h di\ 



+ -x fc x m 9 m a + ^ fei P 2 + -x m d klm P 2 ) 
5 o 9 



1 



r(5) 



1 



— j (4) a vi 



G 



G4)odd = ^ ( i64 3) t/ fe + igHmxt + ^ iW 9 t p 2 ) 



(ffi,)odd 



i£ (3) 

5c 5 « 
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2G 



63 yfe 35 ^ k 



il r 2r(5) 
105 ij 



105 fci 



16 u t(4) 

— x euitiJ.J 
45 R1( -* j)i 



with r 2 = x m x m . The leading components of (g^„) dd 
are identical to those of the 2.5 PN metric in ADM coor- 
dinates. With the help of these formulas, we recover the 
expression proposed by Rezzolla et al. in paper fjj in ab- 
sence of mass multipole moments, i.e. when /^^...i, = 0. 
We conclude that our gauge coincides with the ADM one 
up to the 2.5 PN order, and generalizes that of Rezzolla 
and collaborators at the next level. 

Once the gravitational field has been evaluated in the 
new coordinate system, we may take advantage of the 
freedom in the choice of the matter variables. Vari- 
ables defined by means of the fluid trajectory in the 
configuration space without any reference to the metric 
are particularly appropriate to post-Newtonian calcula- 
tions, since their Poisson brackets are identical to those 
of Newtonian theory [l^. We first take the baryonic 
density p* representing the number of baryon weighed 
by their individual mass per coordinate volume c? 3 x. It 
reduces to the scalar density p at the Newtonian order: 
P* = V~9'P U ° = P + C(l/ c2 )- The second variable we 
shall use for describing the fluid is the coordinate ve- 
locity v l . It does not depend on the metric as being 
the Eulerian quantity associated to the Lagrangian ve- 
locity c&r M (xq , t) I dt of the particle x' i (xg,t) located at 
point Xq at time Xq/c. We have v z /c — u l /u°, and there- 
fore, v l /c — u l + 0(l/c 2 ). We shall introduce a zeroth 
component v° = c when necessary, so that we can write 
v^/c = u^/u . For macroscopic systems, the set of vari- 
ables {p*,v 1 } is completed by the entropy s. 

The relation between a and Ci on one hand, p», v l , 
and s on the other hand, is determined by the form of 
the stress-energy tensor. We shall focus henceforth on 
adiabatic fluid systems, for which 



T* v = p{c 2 + h)u^ + pg 



(4.5) 



where h and p stand for the enthalpy per unit mass and 
for the pressure respectively. The stress-energy tensor 
(|4.5|l has to be supplemented by the equation of state 
of the fluid; it is typically provided by the energy e per 
unit mass 1 as a function of the scalar density and of the 
entropy, e = e(p, s), from what we can infer the pressure 
p(p, s) — p 2 de(p, s) I dp or the enthalpy h(p, s) — e(p, s) + 
p(p, s)/ p. It is now possible to relate the densities a to p*, 
v l , and s. At the Newtonian approximation, we have a = 
(T 00 + T»)/c 2 = [pc 2 u° 2 + O(l/c )]/c 2 = p, + 0(1/ c 2 ), 



1 Let us point out that the symbol e denotes the energy per volume 
unit in paper |^ 



and Oi = T 0l /c = [pc 2 u l u a + 0(l/c)]/c = p*v l + 0(l/c 2 ). 
It is thus useful to pose: 





= G J 


f d 3 y 


P* > 


(4.6a) 






'|x-y| 




=G J 


f d 3 y 


P*v\ 


(4.6b) 






V-y| 




X* 


=G J 


[d 3 y\x- 


y| p* ■ 


(4.6c) 



At the order 1 PN and beyond we proceed in several 
steps, for the computations are longer. We first expand 
\J—g' in power of 1/c in order to evaluate p(p ie ,v p , s) = 
p*/(u°y/—g'). From the expression of the determinant 
of a perturbed Minkowski metric g' — n pl , + Sg' at 
linear order, g' = -1 + Sg' + 0((8g') 2 ) with Sg' = 
&et(ri a p)ri p ' 7 8g' p(J = -if^Sg'^ we get immediately 

^g 1 = (l + ^Sg'^) 1 ' 2 + 0((Sg'f) 

= 1 - ^(<?oo)(<7) + \{g'u){<7) +o(j^\0 (g' odd ) 



h...) + h...) + o 



i 



The lowest contribution to O(l/c 2 )0(g' odd ) is a Lorentz 
scalar proportional to (g' m )( 2 ) or (sy)(2)) and (g' ij )(5)- 
Now, the only allowed combinations are (g la ^)(2) x 
(ffa/u)(5) « 5ij 9ij( 5 y and Tl aP {9' a p)(2)ri p ' 7 {9'pa)(5), which 
are both zero due to the trace-free property of (<?y ) ( 5 j = 
— AG iff /(5c?) . As a consequence, 



S' = 1 + ^ U - 2^o)(7) + 2^)(7) + ? G • •) 

+ 4(-)+ofiV (4.7) 



We next calculate the Lorentz factor vP = dt/dr appear- 
ing in p and if, from the definition of the proper time 
cdr = (-g'apdxO'dxP) 1 / 2 . We find: 



-9oo-29'oi — 



H,3 



-1/2 



i U - (5oo)(7) - 2(ffoi)(6) 



-(g'i 



1 



u 
72 



2 c 2 



l -(...) + h...) + 



-1/2 



(fl r Oo)(7) 



1 v l v^ 
2^0(5) — 



+ i(-) + i(.-) + ofi 



(4.8) 



where i> denotes the Euclidean norm of «*, i.e. w 2 = v % v l . 
Regarding the fact that the lowest odd order part of 
the matter or thermodynamical quantities involved in 
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the stress-energy tensor (|4.5|l is necessarily larger than 
five, equations (|4.7|) and l|4.8l) show that the first odd 
terms in power of 1/c entering the 1/c expansion of 
u° = u°(p*,v p , s) and p — p{p Sf ,v p ,s) arise exactly at 
order seven, which implies that the 1/c 5 parts of the 
scalar density and the Lorentz factor are both zero. Re- 
ferring to these two parts as u ^ roac and p( 5 ) rcac re- 
spectively, so as to distinguish them from the fifth or- 
der of p and u l considered as functions of er, o~i, and 
(Tij, we can write in short: W( 5)rcac = P(5)reac = 0. 
Therefore, the densities a = p(u°) 2 + 0(l/c 2 ) and o~i = 
p{vP) 2 v l + 0(l/c 2 ) have no contribution of order 1/c 5 , 
i.e. <T(5) roac = and <7;(5) rcac = (with the same nota- 
tion as before). Let us mention that the stress density 
0-y = p(u°) 2 v l v j +pg nj +0(l/c 2 ) does contribute at this 
level: 



0"(5)ijroac = P(p* , s) )(5) re 

= -p(p*,s)(g'ij)(5)i 



7^0. 



The fact that we have U, 



(5) rcac 



X(5). 



(4.9) 



and 



i(5) reac 



— for the potentials leads us to conclude that 



the reaction part of the gravitational field is simply 



(ffoo)(7)] 



- 1 



2 I 



+ ^(.--) + G?00)(7) 

§(-Q$*»hu. 



(7) 
rcac 



- [(5oo)(7)] (7) 



G 



-^1 Q$x k d lP < 



(4.10a) 



(#0i)(6) rcac 



+ -{...) + (g' Qi ) m 



= 0. 



(6) rcac 



(4.10b) 



(ffii)(5)] 



^(l + ^)+^(-) + (^)(B) 
[(Sij)(5)](5)rcac 



(5) 



(4.10c) 



at the leading approximation. The functions 
[(ffoo)(7)](7)rcac and [(ff ^ ) (5) ] (5) rcac, are obtained by 
replacing consistently the Newtonian-like potentials 
U, Ui, x, and the mass quadrupole moment 7^, by 
U*, f/*,, X*i and Qij — Jd 3 y y lJ p* respectively in the 
corresponding odd metric components. From equations 
(|4.7|) . JOJ, and (|4. 1U|) . we deduce u° and p expressed 
within our set of variables: 



1 v 2 
2^2 



00 J (7) J (7) reac 



(4.11) 



I I 3 lv 2 



,(s4)(7)](7)reac + [(g'ij ) (5) ] (5) 



reac 9 



i(-)+4(...)+0(i 



(4.12) 



At last, we compute the pressure p as a function of p*, 
v\ and s, by expanding p(p r +Sp, s) = (p* +5p) 2 de(p* + 
Sp 7 s)/dp around p = p*. According to relation JHEJ, 
the difference £p = p — p* splits into (i) a "conserva- 
tive" term of order 1/c 2 , which we shall note P(2) C ons 
henceforth, (ii) a "reaction" term of order 1/c 7 , plus 
(in) contributions of order 1/c 4 , 1/c 6 , and C(l/c 8 ) that 
are irrelevant in our calculation. Thus, C([<5p] 2 ) = 
l/c 4 (. . .) + l/c 6 (. . .) + 0(l/c 8 ), and the quadratic term 
of the Taylor series can be omitted. 

dp 

p(p* + Sp, s) = p(p*, s) + 5p—(p*, s) + O ([Sp] 2 ) 



= Pip*, S) + (p (2 ) + P( 7 )) ^(p*, s) 



13) 



We define the coordinate pressure p* to be the pressure 
of a fictitious fluid of scalar density p* and of entropy 
s in our coordinate grid: p* = p(p„,s). Its numerical 
value a priori differs from that of the actual pressure 
p = p{p{p Sf ,v p ,s),s). We introduce in the same way a 
coordinate enthalpy as well as a coordinate adiabatic 
index 7* = (91np*/<91np*, etc. Whereas h =^ h* and 7 7^ 
7* in general, 7* identifies with the usual adiabatic index 
7 = d lnp/<9 In p when the equation of state is assumed to 
be polytropic. Though 7* happens to be constant in the 
latter case, it usually depends on time and on the field 
point: 7* = 7*(x, t). Within the preceding notation, 
equation l|4.13[) becomes: 



p =p* 



1 - I 3J7* + \v< 



[(flii)(7)](7) rcac + ) (5) ] (5) rcac 



(4.14) 



This entails notably that h = K + (• • -)/c 2 + (. . .)/c 4 + 
C(l/c 6 ). Finally, the 1 PN + 3.5 PN mass density a = 
(T 00 + T M )/c 2 is derived by inserting expressions (|4.3f> . 
(I4.11|) . H4.12p . and (|4.14|) into the stress-energy tensor 
T» v = p*tt°(l + hj 'c 2 )v^v v I '^g 1 ' + pg'^ '. We arrive at 



P* + 



1 



P*C 2 ^1 + U° 2 ) CQns + U° 7) rcac ^ 
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X (l - (v^5 7 )(2)cons ~ (\/-<70(7)reac) ~ P* 
+ ^[P*V 2 +P* (^-[(5'")(5)](5 )r cac)] 

P* + ^2 [ ~ U * + h * + 2 ~ J +/>*( I(5oo)(7)] (7) 
C \ Z P*/ \ reac 



2 [(fit) (7)] (7) reac + 7J )(5)] (5) i 



(4.15) 



noticing that [(# m )( S )](5)reac = -[(Sii)(5)](5)r eac + 
C(l/c 7 ) = C(l/c 7 ). The Poisson integral of the above 
expansion gives the potential U in the {/?*,?/, s} repre- 
sentation, modulo an unimportant prefactor G. 



G 



c^J |x-y 

d 3 y 



x-y 



1 , 



P*[ [(foo)(7)] (7) - q[(5«)(7)] (7) 



+ 2[(.9y)(5)](5)rcac — 



(4.16) 



All other potentials can be cast into a similar form. 
The leading orders of the conservative and reaction part 



are determined by that of a and <x;. As the fifth or- 
der corrections of both mass and current density van- 
ish, we have \ = X* + X(2)cons + X(7)reac and U l = 
U*i + C/j(2) con s + Ui(7)reao modulo irrelevant terms. 

We are now in measure to isolate the "reaction" field 
(ff^)reac- It essentially consists of the odd metric <|4.3|) 
to which add the reaction terms generated by the post- 
Newtonian part of (5^„) C vcn = g'^u - (s^)odd regarded as 
a function of the matter variables a and cr,. The latter 
terms do not contribute to (g' oi ) Teac or (g^reac at the 3.5 
PN order, but enter the 00 component of (gL u )ieac, 

(Soo)rcac = (5oo)odd + ^ U {7) reac + O (-^j > (4.17a) 
(30i)reac = (ffoi)odd - ~j}^i(7) reac + O (jj^J 

= (9oiUd + o(±) , (4T7b) 
2 / 1 

(flij-)reac = (j9ij) odd + 5 %3 ~^ U (7) reac + O I 



(<4) 



odd 







(4.17c) 



By virtue of equations l|4.3jl . I|4.4jl . I|4.16jl . and Ij4.17jl . 

we finally obtain the 2.5 and 3.5 PN components of the 
metric as 



J 



(.9oo)(9) 



(#0i)(8)i 



(fij)(7). 



-A7 



G 



rf 3 y 
x - y| 



d 3 y' 



|y-y 

8G / 1 



+ G 



d 3 y 
x-y| 



?g ) y fc ft(p*« i )-P.« fc QS ) -P*tf fc Q, ( f 



4G 
'5c 7 



f - /, (3) - 20 (3) f/ + — (5) a; fc + -7; fe a^0 (5) - — (5) r 2 - -eu,t 



x 



?3[ -± x k x inM_nW„k 
5 V 21 



(4.18a) 



(4.18b) 



(4.18c) 



with = 



and Sij — Jij(o) cons 
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3.5 PN reaction force 

The evolution equations including the 3.5 PN gravita- 
tional damping are deduced from the conservation of the 
stress-energy tensor 

d a (V=FT£) = IV^FT^d^ , (4.19) 

and from that of the entropy 

d t s + v%s = . (4.20) 

The projection of formula (|4.19|) on it M , together with 
formula l|4.2U[l . leads to the continuity equation 

ftp* + ft(p*<) = 0, (4.21) 

whereas its space part is equivalent to the Euler equation 



P* 



1 + — ) UiV J C 



ft 



-g'T^d t g' a0 . (4.22) 

In the absence of coupling between radiative and non- 
radiative effects, the 1 PN + 3.5 PN dynamics of the 
adiabatic fluid subjected to gravitational damping is de- 
scribed by the post-Newtonian expansion of equalities 
(|4~2UI) . !|OT|) . P2| . including the Newtonian, 1 PN, 
2.5 PN, as well as the 3.5 PN order for fixed p*, v l , and 
s. 

We first identify the momentum density of the fluid 2 
Mi = p»(l + h/c 2 )uic = ^/^g 1 T° i /c on the left-hand side 
of (|4.22(l . The space part Ui of the covariant velocity is 
related to the variables p*, v l , and s through the equality 
Ul = uVc+[«o(9{o +9^/0)] (s+6+8) + (• ■ -)/c 5 + (- ■ -)/c 7 + 
0(l/c 9 ). After replacing m° by its value (|4.11|) in the 
definition of M, and expanding up to the 1 PN + 3.5 PN 
order with the help of equations (|4.3() . we arrive at 



=p*v" + — I h c' - •'>? ", /•' 

P*^ (5y)(5+7) rcac 



.A* 



+ "jP* [h* + -V 2 + U*j V° (flfy)(5) rcac 

+ P*c(goi)(8)rcac 
Note that the momentum density M; is referred to as 7Tj in paper 

s 



P*^ (3oo)(7)roac + (fl'«)(5) 1 



h 1 y 
V V 



h...) + h... ) + o 



(4.23) 



This expression contains the term P*{g'^j)(s,+^)V :, involv- 
ing in particular the third order time derivative of the 
1 PN quadrupole Iij and the fifth time derivative of 
the Newtonian mass multipole moments Qij and Qijk- 
Therefore, the extra time derivation appearing on the 
left-hand side of the Euler equation i|4.22[l induces terms 
in IkL cons , Qf), and Q^f k , which make us lose all benefit 
from the gauge transformation. In order to cure that, we 
modify our system of variables according to paper [f| . We 
keep the baryonic density as well as the entropy, but sub- 
stitute to the coordinate velocity v l the more appropriate 
"momentum velocity" w; = Mi /p* . The relation between 
Wi and v l is found by inverting equation 14.23fl iteratively. 
At the leading order, we have M, = p*v l + 0(l/c 2 ), 
hence Wi — v l + 0(l/c 2 ). This result is used to com- 
pute the 1 PN correction of v l as a function of Wi'. Mi = 
p*v l - p* (—h*Wi - W^w l + A M - w 2 w l /2)/c 2 + 0( l/c 4 ) , 
hence v l , which is in its turn inserted into equation (|4.23|) 
truncated at the next approximation, etc. We get in the 
end: 



v % — Wi H — - ( —h*Wi — 3U*Wi + A»j — -w 2 Wi 



W j(ffy )(5)rcac ~ (ffy ) (7) i 



+ ( 5*7* + ^W 2 + K ) wMj) (5) " ^(9'ij) (5) 

cons rcac 

W k W[ 



XV ' ( 

- C(ff0i) (8) + Y ( -(ffoo) (7) + (9kl) (5) 
rcac \ reac rcat 

+ i(--) + i(--)+0fi) • 



(4.24) 



Up to now, we have reduced the kinetic part of the 
Euler equation l|4.22[) to dj(p r WiV : ') + d t (p*Wi), or equiv- 
alently p^(v J djWi + ftiUj), after the continuity equa- 
tion. The pressure force FP ress = _ di(y/— g'p)/p*, as 
for it, is derived from the post-Newtonian expansion 
(|4.14() and (|4.7|l of p and y '—g' respectively, supple- 
mented by the two approximate equalities U — U* + 
(. . .)/c 2 + (. . .)/c 4 + C(l/c 6 ), (g' u ) {7) = (^)( 7 )reac, and 
v % = Wi — Wj(^j)(5)reac+ irrelevant terms. 



^prcss 



- — di\p* + \ 
P* I c 



2Z7, - 7* ( 3E/» + iw 1 



(d'jj) (7) - (Soo)(7)rcac ~ 7* ( (9jj) (7) + Wjk) (5) 
reac \ rcac rcac 



v 3 v k 



1 / x 1 / 
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2£/* - 



7* (3L7* + IwA) + - ((g'jj) (7) - (g' 00 ) (7) ) - y ( (4) (7) - (&) (5) 

V / / \ rcac reac / \ reac rca< 



WjW k 



dip* 



P* a 
— ft 

P* 

7* 
2 



-=■ 2t/* - 7* 3U* + -w 



((5jj)(7)rcac ~ Ojfe)(5) 



reac 9 



75 ( (Sjj) (7) rcac - (<7oo) (7) ) 



(4.25) 



While both pressure and density tend toward zero in the region of space-time extending outside the system, the ratio 
dip*/p* remains finite by virtue of the thermodynamical relation dip*/p* = dih* — T^diS, with T* = <9e(p„, s) / 'ds 
representing the coordinate temperature. The last term of the Euler equation still to be evaluated is the gravitational 
force i^ grav = dig' a p\/—g'T a P / (2p»). The components g' of the metric in the new gauge and the square root of 
minus the corresponding determinant are given by equations (|4.18|) and l|4.7|l respectively. The stress-energy tensor 
is obtained from the mass density (|4.15|) . the current density a t — p*v % + 0(l/c 2 ), as well as the stress density (|4.9|l . 
by inverting the identities a = (T 00 + T u )/c 2 , a, = T oi /c, and = T ij . All contributions of order 1/c 4 , 1/c 6 , and 
C(l/c 8 ) are neglected as usual. We find: 



-i + - s u*-- i ut 

&■ C 



6 ° I 1 + £ U * 



\J —g' {c 2 <7 — a 

-g'vjk 







c 3 _ 








(-2+0) + 9i 

cons 




V-g'caj 



A* 



-g'caj 



(4) 



(0) 
cons 



2 



-g'(c 2 (T - djj) 



(0) 
cons 

(5) 
rcac 



2 



-g'o jk 



(5) 
rcac 



+ 75^(5oo) (7+9) V-g'ic 2 ' 7 ~ (-2+0) + d i(90j) (6+8) 



-g'ccrj 



(-!)< 



+ 2^(9^) (5+7) 



-9<?jk 



(0+2) c 



o 



1 + y -u*+k+ 2^ 



+ T^C di(p 00 )( 7 ) reac + -OiU* ( (5oo)(7)rcac + (fl'jfc)(5)l 



c 
"2 



^(ff00)(9)rcac + CW J (ff oj ) ( 8 ) 

+ \ \v i v h d i {g' jk ) {7)le&c + ^(4) (7)rcac ) + I(. . .) + I(. . .) + O ( 1 



(4.26) 



The indetermination of the ratio p*/p„ for p» and 
p* both going to zero is raised by use of the re- 
lation p*j p* = p*de{p*, s) / dp*. The replacement 
of v l in Ff r&v according to equation (|4.24|) amounts 
to (i) changing v 1 into Wi at each of its occur- 
rences, and (ii) adding the 3.5 PN term: 5Ff rav — 
-3diU*w J w k (g' :jk ) {5) rcac / c 2 + d l A„ j w k (g' jk ) (5)l . cac /c 2 to 
the resulting expression. Knowing v % , i^ grav 1 and F, ^ prcss j 
we write the Euler equation as: 

d tWl + v j d jWi = if rav + Fr SS ■ (4.27) 

The evolution equations for the fluid are now estab- 
lished. They are parametrized by the Newtonian po- 
tential [/*, the post-Newtonian potential A*, — 4E/*, + 
didtX*/2, as well as the 2.5 and 3.5 PN reaction 

metric Components (Soo)(7) reac, (sdo)(9) rcac, (ffoiWcac, 



(g'ij) (5) reac, and (fir ^ ) (7) reac- Formulas (|4.1U|) and (|4.18|) 
show that these components depend themselves on [/*, 
A*,-, and on the Poisson integrals coming from the func- 
tional derivation. Among these various quantities, the 
only super potential, %*, appears exclusively through 
its gradient, d % x* = Gfd 3 y p*{x % — ?/ l )/|x — y|, which 
can be put under the form diX* = x l A~ 1 (— AirGp*) — 
A^ 1 (—AitGx 1 p*) . Consequently, all the elementary field 
quantities employed in the present formulation are solu- 
tions of decoupled Poisson equations with compact sup- 
port sources. This simple elliptic structure fits particu- 
larly well to numerical computations. The only compli- 
cation, apart from that due to the mere length of the 
expressions, comes from the double integration arising in 
(#oo)(9) rcac- Indeed, the source terms p*Q k 3 ^y k diU* and 
p*A _1 ( —AttGQ^ y k dip-A present in the second member 
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of equation l|4.18bjl are basically themselves Poisson po- 
tentials. This nonlinear effect originates in the fact that 
the integrand of f7(7) rC ac involves certain pieces of the 
2.5 and 3.5 PN gravitational field, namely (<?oo)(7) reao 
(fl , ii)(7)rcac, and (fl^)(5)reac< It occurs either in the orig- 
inal Burke-Thorne-like grid or in the coordinate system 
defined in section ITTll as long as we confine ourselves to 
the {pi,,Wi,s} representation, but can be cured by re- 
defining the density variable as p = p* + ov 7 ). This re- 
sults in a slight modification of the present formalism. 
The corresponding set of equations is presented in Ap- 
pendix^ 

Our last task in the achievement of our initial goal 
will consist in reducing the number of temporal deriva- 
tives remaining in the metric (|4.18(1 to one. This will 
be done by eliminating dtp*., dtWi, and dts with the 
help of the entropy conservation, of the continuity equa- 
tion, and of the Euler equation at the Newtonian or 
post-Newtonian order. The multiple time derivatives 
we shall have to deal with apply on the multipole mo- 
ments. They are of the type (d/dt) n Jd 3 x p*(x, £)/(x, t), 
with n S {2,3,...}. The function / depends on time 
through the evolution variables p*, w%, and s: /(x, t) — 
f(x,p*(y,t),w p (y,t),s(y,t)). Each differentiation is 
performed by means of the formula: 



d_ 

dt 



y^x P »(x,i)/(x,t) = jd 3 Kd t [p^t)f{*,t)] 

d 3 x [-/(x, t)8i (p* (x, t)v* (x, t)) + p* (x, t)d t f{xi, t)] 



d 3 x p*(x,t)J(x, t) 



(4.28) 



where the last equality is obtained by integration by part. 
We recall that the total derivative df /dt = d t f + v l dtf 
represents the temporal variation of the field / for a 
fixed particle of fluid. It satisfies the chain rule df /dt = 
{dp*/dt)df /dp*+(dwi/dt)df /dwi+(ds/dt)df/ds for any 
/ = /(p*, Wi, s). Taking advantage of the Euler equation 
(|4~T7I) . dwi/dt = i f rav + if rcss , as well as the local con- 
servation of mass l|4.21[) . dp*/dt = —di(p„v l ) + v l dip* — 
—p. t diV l , and entropy (|4.2U|) . ds/dt — 0, we get: 



df 
dt 



1L 

9p* 



df_ 



(4.29) 



This shows that the time derivative 



term 



fd 3 



p»(x,t)/(x, t) 



IS 



itself 



of 
of 



the type Jd 3 x p*(x,t)g(x,t) with g = 
(?(x, p*(y, i), w p (y 1 1), s(y, t)). In this way, the number 
of time derivatives acting on the initial integral can be 
brought down to n — 1, n — 2, etc. After n — 1 iteration, 
we are left with a unique time derivative as required. 
The price to pay is an increase of the expression lengths. 
The evaluation of 1^ is particularly tedious since the 
mass quadrupole iy is needed at the order 1/c 2 . After 
specializing the definition (|3.3|l of Zy to the desired level 



of approximation, we have 
1 f'Z 



1 1 j — d x x 

i d 2 



i 



-w 



14c 2 dt 2 
1 



d 3 x r 2 x 13 p 



U* + h 

20 d 
2h?di 



2 P - 

P* 

<i 3 x x ljk p*w k 



O 



We then perform the time derivation along the procedure 
explained above. The force per mass unit i^ +if is 
limited to the Newtonian order: divt/dt — diU*—dip*/p*. 
We are led to Q 



lij — Eij 



Eij -\- Gij 



with 

= STF /d 3 x p*x i x j 



O 



1 A 



(4.30a) 



F i:j = STF /d 3 x 



11 



-w 2 - U* 

7 P* 



p* 

4 .• fc 



— r WiWj --xx w k Wj 
21 7 



+ — x l x j w 2 



(4.30b) 



(4.30c) 



and 



Ga = STF /<i 3 x ^ 



^r 2 x l d U* - ^-x i x 3 x k d k U Jt 
21 21 



(4.30d) 

The computation of 1^ from equation (|4.30a|l is long, 
but straightforward. We make systematical use of formu- 
las (|4~2^|) and The quantities v i and if rav + if rcss 
are truncated at the 1 PN order in the first source term 
of E^, at the Newtonian order in the rest. 



1 ( 1 



dwi 
~dt 



=8iU* 



rcW 2 



r- 2 \ 2™ 



3U, 
+ h, 



1 



r .2 3 



1 



2t/* - 7* ( + 3C/ * ) ) J- 



+ -A 

2— - U* ) diU* 
p* 



p* 



The total time derivative of the enthalpy 
is deduced from relation (|4.29|) : dh^/dt — 
(dp*/dt)(dh*/dp*)(dp*/dp*) = -j*p*d l w l /p* + (D(l/c 2 ). 
For treating the potentials U m , A ti , etc., the operator 
d/dt is first put under its explicit form, d/dt = dt + v l di, 
and the partial derivative dt is next applied on the source 
S through the integral symbol. The source variation d t S 
is transformed by means of the evolution equations at 
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the Newtonian order. For convenience, we may redefine 
as A _1 (— AnGp^Wi). This yields for instance 



—U, = &.- 1 dt(.-*irGp*)+v i diU. 

= - A~ 1 [-4:TrGdi(p*Wi)] + Wid t U* 

-diU^ + WidiU^ + O' ' 



O 



v J djU*i 



=A [— 47rG(— dj(p*WiWj) + p*diU* — <9 4 p*) 
1 



(2) 

As the potential d t {diX*) occurs in I>-\ its time deriva- 
tive must be eliminated as well. In order not to increase 
the number of Poisson equations, we systematically re- 
place its old expression by the new one. It must be de- 
termined at the 2.5 PN order, since A ifi appears in the 1 
PN part of the gravitational force (|4.26(l . 

d t d iX * = d t [x l A-\-^Gp,) - A-^-^GxV*)] 

= -x i d j A- 1 (-4nGp^) 

+ A^i-iTrGx'djip^)} 

= - X i (8 lk - (5j-fc)(5)roac) 9jU* k 

+ A' 1 [-4ttG (S jk - (.g^)( 5)rC ac) x'djfawk)] 

+ \(...) + \(...) + o n 

c z c 1 

(the last form is obtained by using the Newtonian plus 
2.5 PN approximation for the coordinate velocity H4.24fl : 
v % = Wi — Wj(g' i j)^ rcac ). It is sufficient to determine the 

(2) 

total time derivative of d t diX* contributing to I>- ' at the 
Newtonian order. It reads 

9t(d t diX*) 

=d t [—x i djU*j + A~ 1 (—&irGx i dj{p*Wj))\ 
= - x^jA-^-iirGi-dkip.WjWk) 

+ p*d 3 U* — djp*)] 
+ A^ 1 [-ATrGx l (-d jk (p*WjW k ) 

+ dj{p*djU*) - Ap*)] . 

(2) 

After ruling out all temporal differentiation from I>. , 
our result is finalized by removing, whenever possible, the 
space derivatives of the baryonic density p* and the pres- 
sure p* by integration by part. Remark that the space 
and time derivatives of 7, exactly cancel. All quantities 
S\f and Q§1 = (d/dt) n fd 3 x x ijk p* for n < 4 are rewrit- 
ten at Newtonian level following the preceding method. 

Now, Sij and Qijk, as well as the hexadecapole Qijki — 
/d 3 x x^^p*, can be formally assimilated to the multi- 
poles that parametrize the wave form h\^ A in some "ra- 
diative coordinate" grid X M = (T, X) up to the 3.5 or- 
der, as explained in paper After actually achieving 



the substitution into the wave form, the new expression 
for hlj d (X, T) involves one time derivation only, and is 
thus suitable for numerical computations. In conclusion, 
the present formalism allows one to build an accurate 1 
PN + 3.5 PN code (i) that models the fluid dynamics in- 
cluding the 3.5 PN damping effect, and (ii) that permits 
gravitational wave extraction. All the field equations are 
of Poisson type, and have compact support sources. The 
number of time derivatives is never higher than one. 

Let us list in the end the quantities to be computed and 
the equations to be solved for modeling the evolution of 
an adiabatic fluid at the 1 PN + 3.5 PN approximation. 
For reader's convenience, we follow the presentation of 
Blanchet-Damour-Schafer 5|- The notations are as close 
to theirs as possible, though not identical. 3 It is useful 
to pose: 



U5 — (#Oo) (7) rcac , 



(4.31a) 



y 

^7 = y (5oo)(9)rcac + SC^A" 1 {-AllGp, 



* V (5) reac W i ) i 



U5i — -^j-(50i)(8)rcac , 
C 5 

^5ij y \9ij )(5) rcac i 



= y (Sy)(7)reac ~ S ,3 U 5 



The resulting set of formulas is given below. 
Thcrmodynamical quantities: 



(4.31b) 
(4.31c) 

(4.31d) 

(4.31e) 



de(p„,s) 
ds 



_p(p*,s) <9e(p»,s) 

7T* — = yO* 



7 



p* dp* 
d\np(p*,s) 



<91n p* 

/l* = e(p*, s) + 7T* . 

Primary Poisson equations: 

All* = -4ttG P , , 
AU*i = -AirGp^Wt . 

2.5 PN quantities: 

see Appendix IB1 



(4.31f) 
(4.31g) 

(4.31h) 
(4.31i) 



(4.31j) 
(4.31k) 



)[2] 
<ij 



)[2] 



dQ\ 



dt 



We have the following correspondences: 



our notation 


p* 


e 






notation of paper f5] 


r* 


e/r* 


Ui 


A, 



(4.311) 
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h 2G 6 



(4.31m) 
(4.31n) 

,,- r • ■ -/ • ( 4 - 31 °) 

1 PN quantities with radiative corrections included: 
a = 2*7* -7* (^w 2 +3U, 



d t (Q [ k 3 }x i x k d l p 



(4.31dd) 



AR = -4nGQ%x z d jP * , 
U 5 = ^ (-Qflx'djU. + Il 



U 7 = — 



2G 
5 



h [ k lx k diU* + Q k 6 }x k (-diU 2 + 2U,diU.) 



ikl ■ 



~x l d t U* + A H \ + —Q% mX « x l d m U* 
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11 



42 



42 



+ ^ [2C/ 5 + 7* (-3t/ 5 + w k wih 5kl )} , 

/3 = \i 2 + K + 3E7* , 
3 

S = -W 2 +K+ 2tT* - U* , 

Ci = -d 3 (p*WiWj) + p*diU* - c^p* , 
AZ7 2 = -47rGp*<5, 
AD< = -4ttGC 4 , 
AE Z = -iTrGx'djQ , 

2 



(4.31p) 
(4.31q) 

(4.31r) 

(4.31s) 
(4.31t) 
(4.31u) 
(4.31v) 



S$ n x l x n d k U* - 2U*R + R 2 



(4.31ee) 



U 5i 



2G 



1 



\Q^x k d { A xi - \qf k ] A^ k + Q^x'U* 



x l d t R + R4 



(4.31ff) 



2G 



*7y 



k 
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(4.31gg) 



AG = -4ttG <P " - — /i 5jfe x'djip^Wk) , (4.31w) Velocity and forces: 



1 , = 4 [S ik - -%h Uk \ U. k + hji 



1 2 

v' = Wi + — (-(3wi + A^) ^w k h 5ik 



(4.31x) 



4[7 5i + 2w k h 7ik + Wi (3U 5 - w k wih 5kl ) 



3.5 PN quantities: 



2(p + 2U m )w h h 6 ik + 2A* k 



(4.31hh) 



hf}, Q [ S, Q l S, Qgl sg 1 , see AppendixE 



[2] dhf} 



dt 



A[4] _ ^ 
W « " dt ' 

A[4] = d( ?!jfc 
A[4] _ 



Ai? 2 = -4ttG 



(4.31y) 
(4.31z) 
(4.31aa) 
(4.31bb) 

i 2 [ 2 lx k dip, + Q$x k (p&U, + di(pj)) 



if rcss = - (l + ^) (ftft, - T*d iS ) - ^d ia , 

(4.31H) 

+ i \diU 7 - 4w k diU 5k + w k widih 7kl 

+ (U 5 - 5w k wih 5kl ) diU* + (8 + 2U*) diU 5 

+ 2h 5kl w k d l A, l ) . (4.31jj) 



+ 2>p*W k WlQ k 3 } - J^Q k Lx k X l drnP* 

+ Q [ ki xk {^x l d t p* - 4p*wi 

+ Q[1 xk (-^x l x m d mP , + ^r 2 di Plf - p*x l 



5 A 



Evolution system: 

d t p* = -d l (p*v r ) , 
d t s = -v l d l s , 

d tWl = -tfdjWi + F.r v + if 1 ' " ■ 



(4.31kk) 
(4.3111) 
(4.31mm) 



Gravitational wave form (R = VX l X\ N = X/R): 



^ d (X,T) = ^P^(N) \ lif r ^^,„,, 



1 



[3] 



-e klm S$ n x l x n d k p* - 



Ai? ? = -4ttG 



(4.31cc) 

Q§x k d t (p. Wi ) - p.w k Q l $ - p^ fe Q[t ] 



2G 

+ Yc e ^ kS ^ Nb + jh NaNb ^ kl 

+ ^a H kSlf ac N b N c \ . (4.31nn) 
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APPENDIX A: SET OF FORMULAS IN THE 
{p,w ; ,s} REPRESENTATION 

We give here the set of equations describing the 1 + 
3.5 PN dynamics by means of the matter variables p = 
P* + C(7), Wi, and s (see section llV|) . 

Thermodynamical quantities: 



T = 



de{p,s) 
ds 



~ P(p,s) ~de(p,s) 
7T = ; — = p- 



7 



p dp 
<91np(p, s) 



d In p 
h = e(p, s) + tt . 

Primary Poisson equations: 

AU = -AirGp, 
Mfi = -AwGpWi . 

2.5 PN quantities: 

Q^j , see Appendix [5] with the substitutions p* 
U* -> U, etc.; 



(Al) 
(A2) 

(A3) 
(A4) 

(A5) 
(A6) 



AC t 



AttG S j 



c 



5 h 5jk ) x l dj(pwk) ■ 



A,=-\l6 ik -^h 5ik )U k + ~Ci 



1 i I 2 



(A18) 



(A19) 



3.5 PN quantities: 



[2] n [3] n W n [4] q [3] 
Hj ' ^>ij ' ^ij ' ^ijki ij ' 

substitutions p* — > p, t/* — > U, etc.; 



J 2 g J , Qf>, Q™, Q^ k , S\f, see Appendix E with the 



■ [2] = dMj 
2y ' dt 

A [4] _ d Qfj 

w « - dt ' 

dO [4] 

A[4] _ a ^ijk 
^ijk 



s 



[4] 



rfS 



eft 

[4] 



f/f 



2G r 

5 



)[3] 



(x l djU — WiWj) — R 



(A20) 

(A21) 

(A22) 

(A23) 
(A24) 



AR 2 = -4ttG 



i$x k d l p + Q%x*d l (pS) 



+ Apw k wiQ [ k } - j^Q l k Lx k x l d m p 
+ Qlh h (^x l d t ~p-A~pw l 



/ 3] - k - 

5 oM „kj,i 



A [2] d Qfj 

, _ 2G . [2] 


(A7) 
(A8) 


AR t 


Mt=-4nrGQ^x i d j p, 


(A9) 




U 5 = ^(-Qf)x%U + R) . 


(A10) 


u 7 = 


1 PN quantities with radiative corrections included: 




a = 2U - 7 Q™ 2 + 3tA + ^U 5 (l - 7) , 


(All) 




/3 = iw 2 + ft + 3U , 


(A12) 




S = ^w 2 + h + 2n — U , 


(A13) 




Q = -dj(pWiWj) + pdiU - dip , 


(A14) 




AU 2 = -4nGpS , 


(A15) 




ADi = -4ttGQ , 


(A16) 




AEi = -AirGx^Q , 


(A17) 





-e Hm S^ n x l x n d k p 



7, J 



(A25) 



Q [ k 5 jx k d l (pwi) - pw k Qf k ] - px k Qf k ] 



\d t {Q ] §x i x k d l p) 



(A26) 



2G_ 
~5~ 



i 2 [ k }x k d l U + Q k z }x k (-diU 2 



2UdiU\ 



x l d t U + At 



126 

~ \e klm S^ n x l x n d k U - 2UR + R 2 
2G 



-Qty^A.. - -Q%A k + Q%x"U 



>P] A. 



1 



+ -x l d t R + Ri 



(A27) 



(A28) 



20 



tl7ij 



2G 



21 ™ - 2Q™U + ^Qgx k + &x«Q% = ( 1 + 3 J W + + i*^ 



- i*«<«*$** - £*wqn 



21 



IK 



Baryonic density, velocity and forces: 



p* = p ( 1 - j 



(A29) 



(A30) 



If - \ 2 

+ ^ + A i) - -^w k h 5ik 



AU 5i + 2w k h 7lk + Wi (3U 5 - w k wih 5kl ) 



2((3 + 2U)w k h 5ik + 2A k h 5ik ), 



(A31) 



^ipress 



1 + - + - 

c' 



djh - TdiS ] - -^TtdiCt , 



(A32) 



+ ^ ^J7 7 - 4w k diU 5k + w k wid t h 7kl 
(U 5 - 5w kWl h 5kl ) + (s + 2U) diU 5 



+ 



+ 2h 5kl w k d l Aij . 



Evolution system: 

d t p* = -d l (p*v l ) , 
d t s = -v l d l s , 

dtWi = - v j djWi+F r v +Fr ss . 



(A33) 



(A34) 
(A35) 
(A36) 



APPENDIX B: TIME DERIVATIVES OF RELEVANT MULTIPOLE MOMENTS 



This Appendix presents the full explicit expressions for the time derivatives of the multipole moments used in the 
present formalism. 



r(i) _ n [i] , 1 r [i] , n ( 1 

hj -Qij + ^2ij + U [^A 



=STF / <j 2p*x t w ] + - 



(Bla) 



11 2 nrr , 10 „ ^ , 22 2 

- ■ ■ ■ —/>./• 



— X^WiWjWpj 



p*(2A H x j + —r^WjdiU^ + — x p 
+ p*z' [^Y^Wp^pf/* - ^r^djpt/^ - &U*Wj - -^wjw 2 ] 

+ x % x p [p*{-^Wjd p U* + l^WpdjU*) - ^p*(d p Wj + djw p )] + x l X [p*(d p U* p - ^w p d p U*) + ^p*d p w p ] 



17 

+ —p*X l X J X P [d pq U* q - WqdpqU* 



d 3 x + ( - 



(Bib) 



4? = Q l 3 + 4 7 4j + ° ( i ) = STF H W W + 2p*WiWj 



+ 



p*(8D iX j + E iX j + yiZdiUJjU* + AA 



44 
21' 



r 2 Wjd lp U stp — UU^WiWj 
22 , 



+ ^r 2 WjW p d ip U* - ^WiWjw 2 ) - ^r 2 diW p d p Wj + ^(1 - ^)r 2 diWjd p w p + yr 2 %f7*) 



22 2 22 
— 2h*p„,WiWj + 2p*WiW J + —r dih*djp* — —i 



2 T*disdjp* 



12 32 8 8 64 

+ - —WiWjdpU* + —WjW p diU*) + p*(—Wid p Wj + —WjdiW p + -^w p d l w j )] 

23 80 

+ x 1 [p*(2djU2 - 10U*djU* + 10wjdpU* p — 2w p djA* p + 2w p d p A*j + —w 2 djU* — — ttfjiOpdpI/*) 



+ p*(- 2djU* + — 
12 



Wjd p w p 



20 20 
-WpdjW p — TTjWpdpWj) + 2h*p*djU*\ 



21 21 

12 80 32 80 32 

+ x l x p [p*( - —djU*d p U* + 7j-^Wjd pq U* q - -^WpdjqU^q - —WjW q d pq U* + -^w p w q d jq U*) 

4 4 4 4 24 

+ P*(-jd p WqdqWj + -djWqdqWp - -(1 - j*)dpWjd q w q - -(1 - ^)djW p d g w q - —d jp U*) 

8 8 
- -djp^dph* + -T^djPxdps] 



21 



3 47 26 

+ X % X 3 [p,( - -d p U*d p U* + -^WqdpqU^p - —WpWqdpqU*) 

8 8 6 8 8 

+ P*(2y(! - l*){dpWp) 2 - —dqWpdpWq - -(-AnGp*) + —d p Kd p p* - — T*d p sd p p*\ 

11 11 22 11 11 

+ r 2 x % [p*{—d p U*djpU* - -^djpDp - -^Wqdjpqll^p + -^WpWqdjpqU*) + — p*fy(-47rGp*)] 

17 17 34 

+ x z x J x p [p*( - —dqU^dpqU* + -^OpqDq + —w r dp qr U* q 



17 
21 



WqW r 



17 



'*9p(-47rG/9*)] 



rf 3 x + O 



Q'S 



STF y |4p*9jWj + 6p*WjdiU* + 2p. t x l [w p dj p U* - dj p U* p ] j d 3 ? 



[4] 



Q ijk 



Q [1] 

o l3] 



- STF J^Ap^diK - T„dis){djK - T*d jS ) + 2p^(3d i U^d j U„ - A Wj d ip U^ + AwjWpdipU*) 
+ 4p*( — diWpdpWj + diWjdpWp — "f*diWjd p w p + 3dijU* 

+ x l [2/j* (djpU^dpU* - djpDp - 2w q djp q U*p + w p w q djp q U^) + 2p*dj (— 47rGp*)] | d 3 x , 
: STF /|3p*a;Vw fe |rf 3 x, 

= STF J {3p*a;V9 fc t/* + Sp^WjWk } d 3 x, 

: STF y |6p*WjWjtt;fe + a; 1 [l8p*tUfcdjI7* + 12p*djW k ] + 3p*x l x 3 [ - d kp U* p + Wpd kp U*\ j d 3 x , 

: STFfap.vwW. + *8p* Wj 8 iWk + x> [U^djh. - T* dj s)(8 k K - T A s) 

+ p*(l8djU*d k U* - 24w k dj p U* p + 24w k w p djpU*) 
+ p*( — 12djW p dpWk + 12djWkd p Wp — 12^^djW k d p w p + 36dj k U*)] 
+ x l x 3 [p«(3dpU*d kp U* - idkpDp - 6w q d k pqU« p + 3wpW q d k p q U«)] + 3p»dfc(-47rGp*)j d 3 x , 

STF J^p,x i x j x k w^d 3 x, 
STF y \kp*x i x 3 x k d l U if + Uptx'afwkWt} d 3 x , 



Si? = 



,[3] 



+ 

:STF 



(Blc) 
(Bid) 



(Blc) 
(Blf) 

(Big) 
(Blh) 



(Bli) 
(Blj) 
(Blk) 



STF y ^2Ap^x l WjW k wi + x i x j (36p*w l d k U lt + 2Ap*d k w{) + 4p*x i x j x k ( - di p U* p + w p di p U*) j d 3 x , (Bll) 
STF Jep qi {p*x j x p d q U* + p*x p w jWq ^ d 3 x , (Blm) 



STF Je pq i^x p [p^WqdjU* + p*djW q ] 

X P (2p*Wjd q U* + p*d q Wj) + X j p^Wpdqll^ + p*X J X P (w r d qr U* - <9 gr £/* r )J j rf 3 X , (Bill) 

je pqi { [ - p* Wq 8 jWp + p*x p [ - Wq 8 jr U* r + WqWr d jr U*] 

+ p*X P [ — djW r d r W q + djWqd r W r — J+djWqdrWr] ] 

+ [2p*x p (djh* — T^djs){d q h st — T*d q s) + p*w p d q Wj + 3p*WjWpd q U* 

+ 3p*x p [djU it d q U is - Wjd qr U* r + WjW r d qr U*\ + p*x p [d q Wjd r w r - d r w dqW r - "f*d q Wjd r w r + 6dj q U*\ 

+ 2p*X j [wpW r d qr U* - W p d qr U* r \ 

+ x j x p [p*(d qr U*d r U* - d qr D r - 2w s d qrs U* r + w r w s d qrs U*) + p*d q (-AirG p*)\\ | rf 3 x , (Bio) 
STF Jep ql {p lf x J x k x p dqU* + 2p*x j x p w k w q } d 3 x, (Blp) 



22 



Sifk = ^TF j£pqi\2p*x p WjWkW q + p*x j x k w p d q ll* + x j x p [2p*(2w k d q U* + w q d k U*) +2p st (d q w k + d k w q )] 

- p*x 3 x k x p [d qr U* r — w r d qr U*\ | d 3 x . (Blq) 

(Blr) 
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